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The finite size efli'ects on nucleon masses are calculated in relativistic chiral pertur- 
bation theory. Results are compared with two-flavor lattice results. 



1. Introduction 

Finite size effects, in particular on dynamical configurations, can be a se- 
rious impediment to precision lattice calculations of hadron masses and 
matrix elements. On the lattice sizes used in production runs, the nucleons 
and pions are the relevant degrees of freedom for understanding the finite 
size effects on the nucleon mass. We calculate these effects here using two- 
fiavor relativistic baryon chiral perturbation theory (xPT) at 0{q'^) (one 
loop). 

This study of nucleon masses is based on QCDSF and UKQCD un- 
quenched {Nf = 2) lattice data using a Wilson gauge action and two fiavors 
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Figure 1. Nucleon masses from QCDSF and UKQCD (white symbols) and JLQCD 
(gray symbols). 



of non-perturbatively 0(a)-improved Wilson fermions (a denotes the lattice 
spacing). The corresponding pion masses are in the interval 0.45 — 1 GeV. 
Valence and sea quark masses are equal. Lattice sizes are 1 — 2.2 fm. The 
scale is set with tq, using the physical value tq = 0.5 fm ~ l/(395MeV). 
JLQCD has undertaken a study with the same lattice actions and range of 
simulation parameters on varying lattice sizes ^. In the comparison shown 
in Fig. 1 one can observe that masses from lattices smaller than 1.8 fm are 
consistently larger than the results from larger lattices. 

The a ^ limit was not performed. Since discretization effects in un- 
quenched simulations could be quite large ^ we have to keep an eye on 
the discretization errors in our results. Studying the variation of nucleon 
masses with a at a fixed lattice size we find that between two data points 
on approximately the same volume and at the same pion mass, but with a 
varying by ^ 30%, the nucleon mass remains unchanged within the statis- 
tical errors. Moreover, we compare our results from large lattices (L > 1.8 
fm) with two CP-PACS data sets ^ with renormalization group improved 
gauge fields and tree-level tadpole-improved clover quarks, « 1.5 and 
2 GeV respectively and L > 2.5 fm, and find that they agree. 
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Figure 2. Nucleon masses on large lattices. The solid curve represents a fit with rela- 

tivistic xPT, the dashed curve the NR limit using the same values of mg and ci, and 
the dot-dashed curve the NR limit using mo = 0.81 GeV and ci = —1.1 GeV~^. 

2. formalism 

The one-loop contribution at 0{q^) is generated by the 0{q^) Lagrangian 
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U. We use the 



- T Wu^d^u], Ufj, = iu^d^Uu^ and 
infrared regularization scheme ^. To compute the renormalized nucleon 
mass 77? wc include the tree-level 0{q^) term — 4cimp<^, and an additional 
contribution of the form eiWpg, which is needed for renormalization of the 
relativistic 0{q^) result, although formally only entering at 0((/^). The 
one-loop correction to mo is ^: 
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ei{X) denotes the renormalized ei. The renormalization procedure is de- 
tailed in ^ We employ A = IGeV, qa = 1.2 ^ and = 92.4 MeV. 



The nucleon mass in the chiral limit, mo, and the values of ci and are 
determined by a fit to six lattice data points at the smallest masses. The fit 
is shown in Fig. 2. Relativistic xPT describes the data correctly for pion 
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masses up to ~ 700 McV. For the values of the parameters, we find mo = 
0.85(14)GeV, ci = -0.80(18)GeV"^ and e5;(lGeV) = 2.8(1. l)GeV"^ 
These results differ sUghtly from the ones quoted in ^ due to their fitting 
to a larger set of lattice points and using the value gA = 1.267. 

In the limit mps <SC mjv, keeping only terms up to Oinipg) of Eq. (2), 
one obtains non-relativistic (NR) x^^- In Fig- 2, we also plot NR results. 
We find that they approximate the data and the relativistic curve only for 
rather small pion masses. A method to extend the validity of the NR ap- 
proximation to higher momentum scales within a cutoff scheme is described 
in The breakdown of the non-relativistic theory around mps = 400 MeV 
as shown in the plots can be avoided if momentum modes above a cutoff 
Aj^ are absorbed in local counterterms (e.g. see the discussion in ^). Then 
the nonrelativistic 0{q'''} calculation actually agrees quite well with the 
corresponding relativistic result up to mps = 600 MeV. 

3. Finite Size Effects 

We determine the finite size effect from the one-loop 0{q^) contribution to 
the self-energy. Putting the external nucleon line on-shell, the self-energy 
is given by ^ 

S(^ = mo) = 7(2^^^ +"^0^ +mps{i-x)\ (3) 

in Euclidean space. The temporal extent of the lattice is taken to be infinite. 
We define S — ^ (I](^ = too, L) — E(^ = mo, oo)), and using ^ express it 
as an integral over Bessel functions: 

^ = r ^ ^" (L\n\^mlx^+m%s{l-x)) . (4) 

Then we calculate the finite size behavior of the nucleon mass: 

mjv(i') = (1+ <5)mAr(oo). (5) 

We first extrapolate to mjv(oo) by using the lattice result from the largest 
available lattice as input to Eq.(5), and then determine mAr(L) at finite 
lattice extent using the formula Eq. (5) . In Fig. 3 we compare lattice nucleon 
masses from various volumes with both relativistic and NR X-P^ a^t 0{q^). 
Relativistic describes the lattice data very well up to pion masses 

~ 700 MeV. We find that for the finite size effect to be < 1%, L should be 
> 1.9 fm at the smallest pion mass shown in the plot. 
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Figure 3. mjv as a function of size on the lattice and in xPT- Data sets correspond to 
a fixed (/3, k) value and are labeled by the pion mass on the largest volume. The dashed 

lines show the results from NR X-P^- 

The finite size effect calculated in the non-relativistic formalism (for 
discussion see also ^) is considerably smaller compared to the relativistic 
result and the lattice data. On lattices of 1 — 2 £m size it is ~ 30 — 40% of 
the relativistic effect. 
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